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ABSTRACT 

The field of astrobiology has made huge strides in understanding the habitable zones around stars 
(Stellar Habitable Zones) where life can begin, sustain its existence and evolve into complex forms. A 
few studies have extended this idea by modelling galactic-scale habitable zones (Galactic Habitable 
Zones) for our Milky Way and specific elliptical galaxies. However, estimating the habitability for 
galaxies spanning a wide range of physical properties has so far remained an outstanding issue. Here, 
we present a “cosmobiological” framework that allows us to sift through the entire galaxy population 
in the local Universe and answer the question Which type of galaxy is most likely to host complex life 
in the cosmos^^l Interestingly, the three key astrophysical criteria governing habitability (total mass 
in stars, total metal mass and ongoing star formation rate) are found to be intricately linked through 
the “fundamental metallicity relation” as shown by SDSS (Sloan Digital Sky Survey) observations 
of more than a hundred thousand galaxies in the local Universe. Using this relation we show that 
metal-rich, shapeless giant elliptical galaxies at least twice as massive as the Milky Way (with a tenth 
of its star formation rate) can potentially host ten thousand times as many habitable (earth-like) 
planets, making them the most probable “cradles of life” in the Universe. 

Subject headings: Astrobiology, galaxies:fundamental parameters, galaxies:spiral, galaxies:elhptical, 
methods:analytical. Galaxy:fundamental parameters 


1. INTRODUCTION 

A key question in astrobiology is “what is the distri¬ 
bution of habitable conditions and life in the Universe”? 

Observational astronomers have concentrated their ef¬ 
fort on peering into the skies hunting for earth-like “exo¬ 
planets” (see e.g. Seager 2010) and amassed a database 
numbering in the thousand^ In unison, theorists have 
made huge strides in defining the requirements for hab¬ 
itable zones (where planets can produce biospheres) to 
exist around stars of varying ages, masses and metallici- 
tie^ the stellar habitable zones (SHZ; e.g. Huang 1959; 
Hart 1979; Kasting et al. 1993; Rushby et al. 2013; 
Kopparapu et al. 2013; Giiedel et al. 2014). A few theo¬ 
retical works have extended SHZ calculations to estimate 
habitable zones at galactic scales, the Galactic Habitable 
Zone (GHZ; Gonzalez et al. 2001) for specific galaxies 
like the Milky Way (MW) and nearby ellipticals (Suthar 
& McKay 2012; Garigi & Meneses-Goytia 2013). Ap¬ 
proaches to modelling the GHZ of the MW range from 
using metallicity and star formation distributions (Gon¬ 
zalez et al. 2001) to sophisticated descriptions of the dis¬ 
tribution of terrestrial planets (Lineweaver et al. 2004) 
to modelling individual star systems (Gowanlock et al. 
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2011) to tracing the co-evolution of the GHZ and the 
assembly of the MW (and Andromeda) using cosmolog¬ 
ical simulations (Forgan et al. 2015, submitted to Inter¬ 
national Journal of Astrobiology). However, estimating 
the habitability of a whole galaxy population spanning a 
wide range of physical properties has so far remained an 
outstanding issue. 

In this work, we present the first (and simplest) “cos- 
mobiological” (Bernal 1952; Dick 1996) formalism that 
can sift through the entire population of galaxies in the 
local Universe to answer the question Which type of 
galaxy is most habitable in terms of complex life in the 
Universe?^ 

We draw on the progress made by SHZ and GHZ stud¬ 
ies to isolate the primary physical ingredients required 
to model habitability: planets form in stellar vicinities 
with a metallicity-dependent probability since planet- 
formation requires elements heavier than iron. Further, 
explosions of Type H supernovae (SN henceforth) can ex¬ 
pose surfaces of nearby earth-like planets to high doses 
of cosmic radiation, potentially causing partial or mass 
extinctions or delaying the evolution of complex life. Al¬ 
though SN may not sterilise a planet, planets in prox¬ 
imity to regular SN may be less habitable with respect 
to complex life (Clark et al. 1977). We use this under¬ 
standing to propose three key criteria determining the 
existence of habitable oases at the galactic scale: the 
total number of stars, their metallicity and the proba¬ 
bility of planets existing, unmolested by poisonous SN 
radiation. Naturally, galaxies with the largest number of 
potentially habitable planets would be most likely to be 
“habitable”. 
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2. THE MODEL 

We start by linking the habitability criteria mentioned 
above to observable astrophysical quantities: firstly, the 
total number of stars scales with the total mass in stars 
(stellar mass) in a galaxy. Secondly, we make the rea¬ 
sonable assumption that the metallicity of stars closely 
follows that of the gas from which they formed, with 
the gas-phase metallicity being a direct observable for 
galaxies in the local Universe. Thirdly, we make the sim¬ 
plifying assumption that stars are homogeneously dis¬ 
tributed throughout the galaxy so that the probability 
of planets hosting complex life is inversely proportional 
to the fractional volume irradiated by recent SN: a SN 
exploding in a large galaxy will have far less of an im¬ 
pact on far-flung planets as compared to a small galaxy 
where a few SN could expose all of the planets to high 
doses of cosmic radiation, rendering them less hospitable 
for complex life. The fractional volume irradiated by SN 
depends on the volume affected by a single SN multi¬ 
plied by the SN rate, divided by the total volume of the 
galaxy. While some theoretical estimates suggest only 
SN going off closer than 10 parsecs (pcQ to a planet 
can harm complex life or cause mass extinctions (Ellis & 
Schramm 1995), the effects of high radiation doses and 
climatic disturbances on complex life are only poorly un¬ 
derstood at best (Scalo & Wheeler 2002). We encode 
our uncertainty by defining that only SN going off closer 
than Q pc to a planet can potentially harm/inhibit the 
growth of complex life which yields a total affected vol¬ 
ume of |7r((5pc)^. Type II SN explode when the core 
mass exceeds the Chandrasekhar mass limit of 1.44 so¬ 
lar masses producing about 10^^ erg/s of energy, 90% of 
which is radiated away as cosmic radiation (Shklovskii 
1960; Woosley & Janka 2005) that could be potentially 
harmful for complex life. Given that all SN radiate 
roughly the same amount of energy in harmful radiation, 
Q can be assumed to be the same for all SN which results 
in it being factored out in our final equations; the exact 
value of this unknown parameter then never enters into 
our calculations which focus on defining the habitability 
of galaxies relative to the MW. We note that this is an 
upper limit to the total volume irridiated since SN gener¬ 
ally go off in high-density star-formation regions irradi¬ 
ating the same volume. Further, the SN rate is linked to 
the ongoing star formation rate through the initial mass 
function (IMF) that describes the distribution of mass in 
a freshly formed stellar population. For this work, we as¬ 
sume the commonly-used Salpeter IMF (Salpeter 1955) 
for stars in the range of 1 — 100 solar masses (Mq), with 
stars more massive than SMq exploding instantaneously 
as SN; the total volume of a galaxy roughly scales with 
the total stellar mass (e.g. Moster et al. 2010) and we 
ignore the scatter in this relation for simplicity. 

This yields the three key astrophysical criteria that 
determine the habitability of a galaxy : (a) the total 
number of stars (Wot); (b) the probability of these stars 
hosting terrestrial planets or gas-giants {Pt and Pg re¬ 
spectively); and (c) the fraction of these planets that are 
unaffected by SN radiation and therefore have a higher 
probability of supporting complex life (Pcz)- Each of 
these terms is derived as now explained in what follows. 

^ A parsec is a cosmological distance measure equalling 3 X 10^® 
cm. 



Fig. 1. — The fundamental metallicity relation (FMR) linking 
the key physical properties of total stellar mass, star formation 
rate and gas-phase metallicity (Mannucci et al. 2010; Lara-lopez 
et al. 2010). The regions are colour-coded according to the values 
of the absolute metallicity shown by the colour bar. Sloan Digital 
Sky Survey (SDSS) observations show that more than a hundred 
thousand galaxies in the local Universe lie on a very thin plane 
(with minimal dispersion) in this three dimensional space, showing 
an intricately linked physical origin. Our Milky Way (black dot) is 
a “typically average” galaxy, as shown by its position on the plane. 


The total number of stars (Wot) present in a galaxy 
should be proportional to the total stellar mass (M^) 
built up over its entire assembly history such that Wot oc 
M*. 

Next, we link the probability of the galaxy having 
terrestrial {Pt) or gas-giant planets (W), depending on 
its metallicity where we implicitly assume that gas and 
stellar-phase metallicities are correlated. Buchhave et 
al. (2012) have used Kepler data for 226 exoplanets and 
spectroscopy of their host stars to show that small terres¬ 
trial planets (with radius less than 4 times the earth ra¬ 
dius) might be equally probable around stars with metal¬ 
licities ranging between 0.25 — 2.5 solar metallicity (Zq; 
see also Petigura et al. 2013). Given the small sam¬ 
ple, we assume that the probability of finding terrestrial 
planets can be related to the gas-phase metallicity {Zg) 
as Pt oc Zg^ where a is a free parameter. We then study 
Pt in two limiting scenarios: a = 0 {a = 1) such that 
there is no (a strong) probability dependence of terres¬ 
trial planets on the gas-phase metallicity; while the first 
agrees better with observations (Buchhave et al. 2012; 
Petigura et al. 2013), the latter gives an indication of 
the variation expected if the probability were to depend 
strongly on metallicity. To obtain P^, we use the re¬ 
sults of Wang & Fischer (2015) who have obtained high- 
resolution spectroscopy of about 1040 stars in the Keck, 
Lick and Anglo-Australian Telescope planet search sur¬ 
veys to show that the probability of gas-giants scales with 
the metallicity as Pg oc 10^^^. 

Finally, we need to calculate the fraction of planets 
that receive high doses of cosmic radiation due to nearby 
SN. SN explode after roughly 28.6 million years (Myrs; 
Padovani & Matteucci 1993) and so the volume irradi- 
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Fig. 2. — Upper panel: The fraction of total galaxies that can be identified as well-defined spirals (red line) and blob-like ellipticals (blue 
line) as a function of the stellar mass, as observed by the Galaxy And Mass Assembly (GAMA) survey (Kelvin et al. 2014). While low 
mass galaxies are most often spiral systems their fraction drops with mass to roughly half at IO^^Mq with the most massive galaxies 
( ^ 10^ ^Mq) almost all being ellipticals. Lower panel: The number of earth-like habitable planets in galaxies occupying different regions 
of the star formation rate-stellar mass plane, normalised to the MW assuming no-metallicity dependence (o; = 0 in Eqn. 6). The values 
of the ratio (in log) are shown by the colour bar. Our results clearly show that the number of habitable earth-like planets (compared to 
the MW) increases on the diagonal that traces increasing stellar mass and decreasing star formation rate. Low-mass spirals ( ^ 10^ Mq) of 
any star formation rate provide inhospitable environments for earth-like planets to form and evolve. It is predominantly giant ellipticals 
(masses larger than 2.5 times the MW) with low star formation rates (less than a tenth of the MW) that provide the best environment for 
habitable planets to form. 


ated by a SN due to the ongoing star formation {'ip) is Removing the various constants, 


Virr = ^^^^7r(Qpc)^(28.6Myrs) (1) 

where [53.17 Mq]“^ is the SN rate for the chosen Salpeter 
IMF and SN can sufficiently enhance ionizing radiation 
so as to affect planets within a radial value of Q pc as 
explained above. Further, the total volume occupied 
by the stars is taken to be Vtot = M^/p where we as¬ 
sume a homogeneous stellar-mass distribution. The in¬ 
terstellar media of galaxies is assumed to have a density 
p = 200pcrit where pcrit is the critical density of the Uni¬ 
verse. Then, the factional volume irradiated by recent 
SN is 

_ Virr _ ('0/53.17) X (4/3)7r((3pc)^(28.6Myrs) 

( 2 ) 

The probability of planets suitable for hosting complex 
life {Pci) is then inversely proportional to the volume 
irradiated by recent SN as explained above. 

Therefore, the number of habitable terrestrial planets 
will be some fraction of those that have avoided SN irra¬ 
diation, and can be expressed as 

N 

Nt oc NtotPtPci K (3) 

Jirr 

or, 

jv, «_ 0-739M.Z;M./f, _^ 

(V’/53.17) X (4/3)7r(Qpc)3(28.6Myrs) 


Nt^ 


i’ 


( 5 ) 


Using the same approach, the number of habitable gas- 
giants can be expressed as 


,, m2 10^3 


(6) 


Of course, we don’t necessarily mean that the gas-giants 
themselves could be habitable, but that they may host 
moons/satellites that could have conditions suitable for 
life (Kipping et al. 2009). 


3. LINKING HABITABILITY WITH ASTROPHYSIGAL 
OBSERVATIONS 

In an exciting development, observations with the 
Sloan Digital Sky Survey (SDSS) show that the three 
key physical properties (total stellar mass, ongoing star 
formation rate and gas-phase metallicity) are intricately 
linked with minimal scatter for more than a 140,000 
galaxies in the local Universe through an intrinsic “Fun¬ 
damental Metallicity Relation” (FMR; Mannucci et al. 
2010; Lara-lopez et al. 2010) as shown in Fig. The 
same figure shows that the Milky Way is a typical galaxy, 
lying near the knee of the stellar mass-star formation 
rate-metallicity plane with values of 6 x 

^ = 3 M 0 yr“^ and Zg = Zq. It is now understood that 
this relation emerges as a consequence of the ejection of 
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Fig. 3. — The number of earth-like habitable planets in galaxies 
occupying different regions of the star formation rate-stellar mass 
plane, normalised to the MW. The regions are colour-coded accord¬ 
ing to the values of the ratio (in log) shown by the colour bar. The 
two surfaces in this plot are for two different assumptions linking 
the planet-formation probability with the metallicity: the upper 
surface shows results for no-metallicity dependence (i.e. a = 0) in 
Eqn. 6 while the lower surface shows results for a very strong de¬ 
pendence on metallicity (a = 1) in the same equation. As seen, the 
number of planets decreases by about 10,000 assuming a very-high 
metallicity dependence in low mass galaxies. With their high metal 
enrichment, larger galaxies are rather insensitive to metallicity. 

metal-rich gas and the accretion of metal poor gas over 
the entire lifetime of any galaxy (Dave et al. 2012; Dayal 
et al. 2013; Lilly et al. 2013): the position of a galaxy on 
this relation is therefore a “fingerprint” of its entire life- 
history. Given that the metallicity is a function of the 
stellar mass and star formation rate, the total number 
of habitable terrestrial planets and gas-giants can also 
be expressed solely in terms of these two quantities. We 
then obtain the expected number of habitable planets 
for about 140,000 SDSS galaxies by scaling their stellar 
mass and star formation rate compared to the MW as a 
reference. 

Our approach has a two-fold strength: firstly, as a con¬ 
sequence of using the FMR our model implicitly includes 
the entire formation history of all the galaxies in the local 
Universe, a critical necessity for modelling GHZ (Gonza¬ 
lez et al. 2001; Lineweaver et al. 2004; Gowanlock et al. 
2011). Secondly, by using the MW as a “reference” for 
habitability, we are able to remove all the constants of 
proportionality to obtain Eqns. 6 and 7. 

4. RESULTS AND DISCUSSION 

Our results show that compared to the MW, the num¬ 
ber of habitable terrestrial planets increases with increas¬ 
ing stellar mass and decreasing star formation rate^ as¬ 
suming the simplest case of no-metallicity dependence 
(i.e. a = 0 in Eqn. 1), as shown in Eig. [3 Physically, 
this is sensible because the number of stars (and hence 
planets) increases with the stellar mass, with a low star 
formation ensuring most of these systems are unaffected 
by harmful SN radiation. On the other hand, low mass 
galaxies will host a comparatively small number of plan¬ 
ets, most of which will be rendered less hospitable by 
SN. To quantify, low-mass galaxies ( 10^ solar masses) 

with any star formation rate provide less hospitable envi¬ 
ronments compared to the MW, typically hosting a hun¬ 


dred to ten-million times less earth-like planets capable 
of supporting complex life. Its is predominantly massive 
galaxies (masses larger than twice the MW) with less 
than a tenth of the star formation rates that provide the 
best environment for habitable planets hosting complex 
life. 

We then use results from the Galaxy And Mass As¬ 
sembly (GAMA) survey (Kelvin et al. 2014) which is 
the largest study linking the stellar mass to the mor¬ 
phology (shape) of local Universe galaxies. The results 
of this survey (upper panel of Eig. show that low 
mass ( ^ 10^ solar masses) galaxies are most often spi¬ 
ral systems with a fraction that drops to roughly half at 
a stellar mass of 10^^ solar masses. However, the most 
massive galaxies ( ^ 10^^ Mq) are predominantly all el¬ 
lipticals, with a fraction that drops steeply with decreas¬ 
ing mass. Therefore our results show that giant ellip¬ 
ticals with mass ^ and star formation rates 

^ O.SMq yr~^ offer the most hospitable environments 
for planets to form and can host 100-10000 times more 
earth-like planets than the MW. 

The lowest mass galaxies (M^ ^ 10^ solar masses) with 
the largest star formation rates ( ^ 1 solar mass per year) 
are the most metal poor in the local Universe as a re¬ 
sult of SN pushing out most of their metals. A strong 
metallicity dependence (a = 1 in Eqn. 1) has the largest 
impact on these tiny metal-depleted systems: while the 
the number of habitable terrestrial planets decreases by 
about 10,000 for these lowest mass galaxies as compared 
to no-metallicity dependence, massive metal-rich galax¬ 
ies are effectively insensitive to the metallicity term, as 
shown in Eig. 

Similar to the behaviour shown for earth-like planets, 
low-mass spirals ( ^ 10^ Mq) of any star formation rate 
provide less hospitable environments for gas-giants to 
form and evolve as shown in Eig. It is predominantly 
giant ellipticals (masses larger than twice the MW) with 
a tenth of its star formation rates that provide the best 
environment for gas-giants that can host complex life. 
Indeed, the most-massive lowest-star forming giant ellip¬ 
tical galaxies can host as many as 1,000,000 gas-giants 
compared to the MW. 

To conclude, our quest was to answer the question, 
“Which type of galaxy is most habitable”? Sifting 
through a sample of more than 140,000 local Universe 
galaxies, we find that the most habitable galaxies are 
giant ellipticals with masses larger than twice that, and 
star formation rates less than a tenth of the MW. These 
shapeless cosmic giants can host ten-thousand times as 
many earth-like planets and a million times as many gas- 
giants capable of supporting complex life as compared to 
the MW, making them the sites of the largest number of 
potentially habitable or life-bearing planets in the Uni¬ 
verse. 

The closest giant elliptical to the milky Way is Maffeil 
(Maffei 1968) which was discovered in 1968 by the Italian 
astronomer Paolo Maffei. The latest measurements (Ein- 
gerhut et al. 2003) suggest this galaxy lies at a distance 
of 2.92 Mpc (where a Mpc = 3.086 x lO^^cm) from our 
galaxy: if we were listening to our cosmic neighbours, any 
light signal transmitted from Maffeil would take about 
9.2 Million years to reach us. 

We end by noting that our endeavour has been to 
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Fig. 4. — Upper panel: The fraction of total galaxies that can be identified as well-defined spirals (red line) and blob-like ellipticals (blue 
line) as a function of the stellar mass, as observed by the GAMA survey and explained in Fig. Lower panel: The number of Neptune-like 
(gas giant) planets in galaxies occupying different regions of the star formation rate-stellar mass plane, normalised to the MW. The regions 
are color-coded according to the values of the ratio (in log) shown by the colour bar. As for earth-like planets, low-mass spirals ( ^ 10^ Mq) 
with any star formation rate provide inhospitable environments for gas-giants to form and evolve. It is predominantly giant ellipticals 
(masses larger than 2.5 times the MW) with low star formation rates (less than a tenth of the MW) that provide the best environment for 
gas-giants to form. 


present the simplest cosmobiological formalism to hunt 
for the most habitable type of galaxy. Whilst our ap¬ 
proach accounts for the entire formation history of all the 
galaxies in the local Universe, a critical necessity for mod¬ 
elling GHZ, we have naturally made a number of sim¬ 
plifying astrophysical assumptions and approximations: 
we assume stars (and hence SN) are homogeneously dis¬ 
tributed in galaxies, the gas and stellar metallicities trace 
each-other and that the volume scales with the total stel¬ 
lar mass. Finally, the effects of SN cosmic radiation in 


causing partial or mass extinctions or delaying the evolu¬ 
tion of complex life are only poorly understood at best. 
We recognise all these assumptions have ample scope for 
further improvement which we aim to tackle in future 
studies. 
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